Ab-initio prediction of a new multiferroic with large polarization and magnetization 
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We describe the design of a new magnetic ferroelectric with large spontaneous magnetization and 
polarization using first-principles density functional theory. The usual difficulties associated with the 
production of robustly-insulating ferromagnets are circumvented by incorporating the magnetism 
through /erri-magnetic behavior. We show that the the ordered perovskite Bi2FeCr06 will have a 
polarization of ~80 /iC/cm 2 , a piezoelectric coefficient of 283 fiC/cm 2 , and a magnetization of ~160 
emu/cm (2 fj,s per formula unit), far exceeding the properties of any known multiferroic. 
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Multifunctional materials that combine a spontaneous 
magnetization with a ferroelectric polarization are of 
tremendous fundamental and practical mteves^i^^. For 
example they have high dielectric permittivity and high 
magnetic permeability, and could therefore replace the 
inductor and capacitor in resonant circuits with a sin- 
gle component, further miniaturizing portable cellular 
technologies. Strong coupling between the polarization 
and magnetization would allow ferroelectric data stor- 
age combined with a magnetic read. And the ability to 
tune or switch the magnetic properties with an electric 
field and vice versa could lead to as-yet-unanticipated de- 
velopments in conventional devices such as transducers. 
However a single phase material with large and robust 
magnetization and polarization has not been previously 
identified; known ferromagnetic ferroelectrics tend to 
have low magnetic Curie temperatures^, are often weak 
ferromagnets^, or are not strong enough insulators to sus- 
tain a ferroelectric polarization at room temperature^. 

The coexistence of ferromagnetism and ferroelectric- 
ity is difficult to achieve for many reasons^. First, and 
most fundamentally, the off-centering of transition metal 
ions which creates the electric polarization in conven- 
tional ferroelectrics is driven by a second-order Jahn- 
Teller distortion which requires a formally empty d- 
electron configuration 8 . In contrast, ferromagnetism re- 
quires unpaired electrons, which in many materials are 
provided by d electrons on transition metal ions. So the 
coexistence of the two phenomena, although not prohib- 
ited by any physical law or symmetry consideration, is 
discouraged by the local chemistry that favors one or 
the other but not both. In practice, alternative mecha- 
nisms for introducing off-centering on ions other than the 
magnetic transition metal ions, such as stereochemically- 
active lone pairsi^ or geometrically-driven distortions 11 , 
can be used to circumvent this restriction. Second, and 
important from a practical point of view, is the fact 
that ferroelectrics must be robust insulators, so that an 
applied field can reorient the spontaneous polarization 
without causing conductivity. In practice, many ferro- 



magnetic materials are metallic, and magnetic insula- 
tors tend to show aniiferromagnetic coupling between 
the magnetic ions. Finally, magnetic transition metal 
ions tend to be fairly easy to oxidize or reduce, and are 
associated with multiple valence states accompanied by 
anion non-stoichiometry. So those few ferromagnetic in- 
sulators that do exist in fact show hopping conductivity 
at room temperature^. 

In this work we propose a new class of robustly- 
insulating magnetic ferroelectrics in which all three of 
these limitations are circumvented. We use the 6s 2 lone 
pair on Bi 3+ ions, well-established to be the source of 
the ferroelectricity in the multiferroics BiMnOgi^ and 
BiFeOgi^, to introduce the off-center distortion. To avoid 
the difficulty of finding insulating ferromagnetic behav- 
ior, we instead choose to introduce spontaneous magneti- 
zation via ferrimagnetism. In ferrimagnets, the coupling 
between neighboring magnetic moments is robustly an- 
tiferromagnetic, but two antiferromagnetic sublattices of 
different magnetizations are incorporated to provide a 
net magnetization. Because the antiferromagnetic su- 
perexchange is so strong, the magnetic Curie tempera- 
tures in ferrimagnets are usually far above room temper- 
ature. In addition, to promote strongly insulating behav- 
ior, we choose transition metal ions which are resistant 
to oxidation and reduction: d 3 Cr 3+ , in which the up- 
spin t2 S -shell is filled, and high-spin d 5 Fe 3+ , with its 
completely filled up-spin manifold. 

Having selected our cationic composition (Bi 3+ , Cr 3+ , 
Fe 3+ ), our next step is to choose a trial arrangement of 
the ions, and to calculate the crystal structure and mag- 
netic properties. Motivated by the report of the possi- 
bility of single-atomic-layer superlattice ordering in dou- 
ble perovskite L^FeCrOe^, we choose as our trial com- 
pound the ordered perovskite Bi2FeCrOg, with planes of 
Fe and Cr ions alternating in the [111] direction; this 
corresponds to a rocksalt ordering of the Fe- and Cr- 
octahedra, as shown in Figure (a) . We use density 
functional theory within the LDA+U methodic, as im- 
plemented in the VASP package^, to determine the prop- 
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erties. This extension to the usual local spin density 
approximation (LSDA) is necessary, because the LSDA 
gives qualitatively incorrect results for systems with a 
strong local correlation. In our case, LSDA predicts 
that Bi2FeCrC>6 should be a metal containing low-spin Fc 
ions, instead of the high-spin insulator obtained within 
the LDA+C7 method. Technical details include a U 
value of 4 eV, a J value of 0.8 eV, a 6x6x6 k-point 
grid for the two formula-unit rhombohedral unit cell, the 
projector-augmented plane wave (PAW) method with de- 
fault pseudopotentials 16 and an energy cut-off of 450 eV. 
[7=4 eV is the smallest value that gives an insulating 
centrosymmctric structure; this represents a lower limit 
for U, which is not inconsistent with the experimental 
estimation for LaCrC>3 and LaFeO^Ai. 

First we calculate the lowest-energy crystal structure 
using the standard method of adjusting the positions of 
the atoms to minimize their Hellmann-Feynman forces. 
To ensure that our calculations explore all likely sym- 
metries in a reasonable computational time, we use a 
range of starting structures for the optimization, which 
we obtain by freezing in combinations of the unstable 
phonon modes of BiCrOgi^. We find that the lowest en- 
ergy structure has R3 symmetry (spacegroup 146) and 
exhibits alternating rotations of the oxygen octahedra 
along the [111] direction, combined with relative displace- 
ments of the anions and cations along [111] (see Fig- 
ure ^) . This symmetry permits the existence of ferro- 
electricity, and indeed we calculate a ferroelectric polar- 
ization (using the Berry phase formalismi2i22i2i) of 79.6 
/xC/cm 2 . The ferroelectric distortion is driven by the 
well-established stereochemical activity on the Bi lone 
pair— as can be seen clearly in the calculated electron 
localization function 22 (Figure [2J. We calculate a piezo- 
electric coefficient along [111] of 283 /iC/cm 2 , comparable 
to that of PbTiOj 23 in which the piezoelectric coefficient 
along the polar z axis, 633, is 323 /iC/cm 2 . 

Next we determine the magnetic ordering. We test 
all the likely possibilities that are accessible in our com- 
putations: ferromagnetic and antiferromagnetic ordering 
of the Fe planes and Cr planes, coupled both ferromag- 
netically and antiferromagnetically to each other. (Note 
that our super-cell technique does not allow us to study 
less likely long- wavelength magnetic orderings, and so we 
have not explored more obscure magnetic arrangements, 
or non-collinear spin structures). Here we find that the 
most favorable magnetic ordering consists of ferromag- 
netic planes of Fe and Cr, that are antiferromagnetically 
coupled to each other. Since the magnetic moments of 
Fe and Cr arc different (5 /.ib for Fe 3+ and 3 hb for 
Cr 3+ ) this arrangement results in a net magnetic order- 
ing of 2 hb per Fe-Cr pair, which corresponds to 160.53 
emu/cm 3 . 

Thus we predict that the properties of Bi 2 FeCr06 will 
surpass all previous reports for coexisting magnetism and 
ferroelectricity. The ferroelectric polarization is com- 
parable to the best known ferroelectrics (the value for 
PbTiC>3, for example, is ~ 75 /j,C/cm 2 — ), the piezoelec- 



tric response is substantial, and the large atomic displace- 
ments and polarization reflect a high ferroelectric Curie 
temperature. The material is insulating; we have chosen 
transition metal ionic configurations that are the most 
robust to oxidation and reduction. And the spontaneous 
magnetization is comparable to that of the cubic ferrites 
(for example the magnetization of nickel ferrite is ~ 270 
emu/cm 3 ). Finally we mention that, although we de- 
liberately designed our system to be a robust insulator, 
if carriers are introduced either by doping or by optical 
excitation, the transport and optical properties could be 
useful for spintronic applications. The calculated density 
of states, shown in Figure |3 indicates that both the top 
of the valence band and the bottom of the conduction 
band are 100% spin-polarized, with electrons confined to 
narrow, up-spin Fe 3d bands, and holes in broader, up- 
spin Cr 3d - O 2p hybrids. 

Our work illustrates the viability of electronic struc- 
ture methods in the design of new and multifunctional 
materials and offers a synthetic challenge which we hope 
will stimulate further research on this emerging class of 
multiferroic materials. 
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FIG. 1: Comparison of ideal cubic perovskite-structure 
Bi2FeCrOe (a) with our calculated ground state structure (b), 
looking down the y axis of the ideal cubic unit cell. (Ill) 
planes of Fe ions (light gray) alternate with planes of Cr ions 
(dark gray) along the [111] direction. Each transition metal is 
surrounded by an octahedron of oxygen anions (black), and 
the Bi ions (medium gray) occupy the corners of the ideal 
cubic unit cells. The transition from the cubic to the ground 
state structure consists of two distortions: rigid rotations of 
adjacent oxygen octahedra in opposite directions around the 
[111] axis, and displacements of the anion cages relative to 
the cations in the [111] direction. 




FIG. 2: Calculated electron localization function (ELF), at a 
value of 0.65, for Bi2FeCr06- The polar [111] axis is oriented 
vertically (the Bi - Fe - Bi - Cr - Bi ordering can be seen clearly 
along the front right edge of the cell). The lobes are regions 
of high electron localization associated with the lone pairs 
of electrons on the Bi ions. A contour plot of the electron 
localization in a slice at the back of the cell is also shown. 
The ELF was generated using the STUTTGART TB-LMTO- 
ASA-cod<*25 with the LSDA formalism. 
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FIG. 3: Calculated density of states (DOS) of Bi 2 FeCr0 6 . 
The top of the valence band is set to eV. The black line 
shows the total density of states, and the dark and light shad- 
ing the contributions from the Fe 3d and Cr 3d states respec- 
tively. The states at ~ -10 eV are the Bi 6s states, and the 
unshaded states in the broad band below the Fermi energy 
have largely O 2p character. 



